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Abstract. We present some properties of measures (g— Gaussian) that or- 
thogonalize the set of g— Hermite polynomials. We also present an algorithm 
for simulating i.i.d. sequences of random variables having g— Gaussian distri- 
bution. 



1. Introduction 

The paper is devoted to recollection of known and presentation of some new prop- 
erties of a distribution called Gaussian. We propose also a method of simulation 
of i.i.d. sequences drown from it. 

(7— Gaussian is in fact a family of distributions indexed by a parameter q G [—1, 1]. 
It is defined as follows. 

For g = —1, it is a discrete 2 point distribution, which assigns values 1/2 to —1 and 
1. 

For q Cz (—1,1), it has density given by 



00 00 



2^v/4-(l-9) 



A;=0 k=0 



for |a:| < ■ In particular fn ix\0) = ^-\/4 — x^, for |a;| < 2. Hence it is Wigner 
distribution with radius 2. 

For q = I, q- Gaussian distribution is the Normal distribution with parameters 
and 1. 

Below, we present plots of fn {x\ — .4) in blue, fn ix\.l) in orange, Jh {x\-^) in 
red and standard normal density in black 
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This family of distributions was defined first in the paper of M. Bozejko, B. 
Kiimmerer and R. Speichcr in 1997 in in nonconimutative probabihty context. 
Later ( 8 ) it appeared in quite classical context namely as a stationary distri- 
bution Ph of discrete time random field X={X„}^g2 defined by the following 
relationships: E (X,) = 0, E {Xf) = 1, i e Z, 

(1.1) 3a e R;Vn e Z : E(X„|JP^„) = a(X„_i +X„+i), a.s. 
and 

(1.2) 3yl,B,CeM;VneZ:E(x2|Jc-_^„) = 

(1.3) A {Xl_^ + Xl+^) + BXn-iX^+i + C, a.s., 

where T^m '-^ cr [Xk '■ k ^ m). It turns out that parameters a, A, B, C are related 
to one another in such a way that there are two parameters q > and < |p| < 1 
and all others can be expressed through them: 

l+p2' (p2 + 1) (1 _ gp4) ' 

(p2 + 1) (1 - qp^) ' l-qp^ ■ 

Then, one proves that 

(1.4) Vn e Z,A:,z > 1 : E (ff,. (X„|g) |.F<„_,) = p'^'Hk (X^-M , a.s., 

where J-<m f [Xk ■ k < m), (similarly one defines J->m ■— c (^fe '■ k > m) ) and 
Hk {x\q) : k > ~1 are g— Hermite polynomials defined below. It turns out that for 
q > I the one-dimensional distribution of the process X is not defined by moments. 
This case is treated separately (e.g. in in [12] ). 

As mentioned earlier, here we will consider only the case |g| < 1. We will preserve 
notation and denote family of g— Gaussian distributions by Ph (?) or simply Ph. 

It turns out that there is quite large literature where this distribution appears 
and is used to model different phenomena. See e.g. [H], [5], [2], [6], 
[9j. Besides random field defined above models notions that first appeared in non- 
commutative context and hence establishes a link between noncommutative and 
classical probability theories. 

Remark 1. In the literature there exists another family of distributions under the 
same name. It appears in the context of (Boltzmann-Gibbs)-statistical mechanics. 
See e.^.flS] for applications and review. 
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Both families are indexed by basically one parameter q G [—1, 1], and for q — 1 
both include ordinary N (0, 1) distribution. 

In the sequel we will use the following traditional notation used in so called '9- 

n 

series theory' [0]^ = 0, [n], = 1 + g + . . . + n > 1, [0]^! = 1, [n]^! = Y[ [i]^, 

1=1 

for n > 0, k = 0, ...,n and otherwise, and {a\q)^ = 



k 

rn-l 



[k]Un-k] 



Y[k=Q (1 ~ '^1'') , ioT n — 1,2, . . . ,00 (so called Pochhammer symbol). Sometimes 
i^\l)n ^^^1 abbreviated to (a)^ if it will not cause misunderstanding. Notice that 

(g)„ ^ (1 - qr N,!, C], = TOtTT that , [n]^! and 

! and (^) (Newton's symbol) respectively as > 1. 



n 
k 



tend to n, 

9 



Remark 2. Introducing new variable z defined by the relationship 1 — ^ "^-^^ = 
cos 2712: we can express q^ Gaussian density through Jacobi 9 functions defined e.g. 
in |14j . Namely we have for z G [—1/2, 1/2] 

/H(2sin(7rz)/v/r^|g) =C,03 (219)^2 (z|g) , 

TOt/i Cq = 2')rgi/-t(g2 j^ wkcre 63 {z\q) and 62 {z\q) are so called third and second 
Jacobi Theta functions. 

Let us introduce family of polynomials (called 5— Hermite) satisfying the follow- 
ing three-term recurrence relationship 

(1.5) Hn+i{x\q) = xHn{x\q) - [n]qHn~i{x\q), 

with H_i{x\q) = 0, Hq {x\q) = 1. Notice that i7„ (a;|0) = C/„ (x/2) ,n>~l, where 
Un [x) are Chebyshev polynomials of the second kind defined by 

(1.6) Un {cosO) = sin((n + l)6')/sin6i, 

and that Hn{x\l) — Hn{x) , n > — 1 where Hn{x) are ('probabilist') Hermite 
polynomials i.e. polynomials orthogonal with respect to Gaussian N (0, 1) measure. 

It turns out that q- Gaussian is the distribution with respect to which q-Hermite 
polynomials are orthogonal. This fact can be easily deduced from (|1.4p . 

Thus in particular using the condition 

Vn> l;ge (-1,1) : / {x\q) fn {x\q) dx ^ 0, 

we can get all moments of q- Gaussian distribution. Hence in particular we have 
]g^2n+i ^ 0, n > 0, ^ 1, ^X^ = 2 + q, EX*^ = 5 + 6q + 3q^ + q^ if only 

X^Ph (q) ■ 

The aim of this paper is to make q- Gaussian distribution more friendly by pre- 
senting an alternative form of the density fn for q e (—1,1), more easy to deal 
with (in particular we find the c.d.f. of Ph), and suggest a method of simulation 
of i.i.d. sequences having density fn- 
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2. Expansion of /h 

In this section we will prove the following expansion theorem: 

Theorem 1. For all \q\ < 1, |a;| < — q we have 

/I °° 

(2.1) fHix\q) = ^^^y4^7T^^^(-l)fe-iq(^)[/2fe_2 [x^T~q/2 

k=l 
where (l) = 

As a corollary we get expression for the c.d.f. function of Ph- 
Corollary 1. The distribution function of q— Gaussian distribution is given by 

(2.2a) 



(2.2b) 



Fh {y\q) = \ + ^ arcsin (^yy^l - q/2 



2tt 



k=l 



2k 



Identity (12. ip can be a source of many interesting identities, which may not be 
widely known outside the circle of researchers working in special functions. 



Corollary 2. For all q e (-1, 1) 
-9' VT^J 



(2.3) (9)^ n ((1 + - E(-l)'"'9(')c/2.-2 {x^/T^/2) 



fe=i 



where the polynomials Uk are defined by U.6\) . 
In particular, we have: 

a) 



k=l 



a particular case of so called Jacobi's 'triple product identity'. 
Hi) 



iv) 
(2.4) 



= 1 



k=l 



('?)L = i + E(-i)'+M2fc-i)g(^), 



k=2 



(2.5) 



Vn e N 



2n 
n 



9 k=l 



2n 
n — k 



vi) fH{x\q) is bimodal for q G (—1,(70)7 where (7o(— —.107) is the largest real 
root of the equation E^o(2fc + l)2g'=('^+i)/2 = 0. 
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Lemma 1. For all q G (—1, 1) and n> A we have 
i) 



sup 



< !M1 



fH {x\q) ~ V4 - (1 - 9) x2(^(-l)'=-lgG)c/2fc_2 [x^T~q/2 

(n-l)(n-2)/2 



fc=l 



sup \FH{y\q)- t:. arcsin ( - q/2 



fc=i 

|^|n(n-l)/2 ^ |^|n(n+l)/2 

Remark 3. Using the assertion of the above corollary one can approximate the 
density fn as well as function Fh {v\q) ~ 5 ~ :^ arcsin {ijy/l — q/2) of Ph by ex- 
pressions of the type -^4 — (1 — g) x polynomial in x with great accuracy. This 
expression is simple to analyze, simulate and calculate interesting characteristics. 
Of course one should be aware that for small values of n ^^^"^ -\/4 — (1 — g) x 

(ELi(^l)''""^9*''''^2fc-2 is not nonnegative for all (1 - q)x^ < 4/ To 

give a scent of how many n 's are needed to obtain the given accuracy we solved 
numerically (using program Mathematica) the equation 

I |(ri-l)(n-2)/2 



7r(l-g2)2 

for several values q and e. Let us denote by N [q, e) the solution of this equation. 

e=\q= .1 .4 .7 .9 .99. 
, .01 3.59 4.97 7.71 14.93 56.73 

.001 4.04 5.67 8.73 16.53 60.86 ' 

.0001 4.76 6.26 9.61 17.97 64.70 
We also performed similar calculations for equation 

, |n(n-l)/2 _^ , |n(n+l)/2 



obtaining: 









27r(l 




) 


e ^\q = 


.1 


.4 


.7 


.9 


.99 


.01 


2.3 


3.3 


5.1 


9.5 


33 


.001 


2.8 


4.2 


6.3 


11.5 


39 


.0001 


3.2 


4.73 


7.3 


13.3 


44 



3. Simulation 

There is an interesting problem of quick simulation of i.i.d. sequences drawn 
from q — Gaussian distribution, using few realizations of i.i.d. standard uniform 
variates. One possibility is the rejection method (see for example [IH])- It is not 
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optimal in the sense that it uses least realizations of independent, uniform on [0, 1] 
variates. But as one can see below it works. 

To apply this method one has to compare density of the generated variates 
with another density that has the property of being 'easy generated' or another 
words i.i.d. sequences of variables having this control density are easily obtainable. 

/ {\ — q\U — (\ — q\x^) 

In the case of density Jh such natural candidate is — . However 

this density is unimodal, while the densities Jh for q below certain negative value 
are bimodal. This would lead to inefhcient simulation method requiring many 

trial observations to be generated from ''^^ ■ to obtain one observation 

from fn for sufficiently small q. That is why we decided to take as 'easy' density 
the following one: 

V(i - g)(4 - (1 - q)x^) n,Li((i + q'f - (1 - qW^^) 

^"'^^ ^ Mm ' 

defined for a; e ( ,7^ , —i — ) . However to be sure that this distribution can be 
used one has to prove the following inequalities presented by the following Lemma. 



Lemma 2. For — 1 < g < 1 and x E 



1-9' VT^ 

Ih {x\q) 



we have: 



Ie {x\q) 

where 

( (l + g)(l-g3)(l-g'^)(l-99)nr=4(l-'?'')(l + '?') 9 e (0, 1) 

^^^'^"1 (i+9)(i-9-')(i-9')(i-9^)nr=4(i + M')'(i-'?') ^./'zec-i^o) 

Function M (q) has the following plot 



Fig. 1. M{q){0 ^ 5) versus q (-1 1). 

Now following [lOj we can simulate sequences of independent random variables 
with q — Gaussian distribution. If g = ±1 then such simulation is trivial. 

For q E (—1,1) we use Lemma [2] and program Mathematica. We generated 
sequence of independent random variables from density fs by inversion method (see 
[lOj). since fE{x\q) can be integrated leading to cumulative distribution function 
(c.d.f.): 

/ fE iy\q) dy^D (q) x^{l -q){A-{l- q)x^)Q2 {x\ q) 

+ C {q) arcsin(xv/r^/2) Fe {x\q) , 
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for I a; I < ^i-q i where Q denotes quadratic polynomial in x'^ with coefficients 
depending on g, while the constants D and C are known functions of q. Recall 
that the inversion method requires solving numerically the sequence of equations 
Fe {x\q) = Ti, where n are observations drown from standard uniform distribution. 

Since the function Fe {x\q) is strictly increasing on its support and its derivative 
is known, there are no numerical problems in solving this equation. Due to efficient 
procedure 'FindRoot' of Mathematica solving this equation is quick. 

Now let us recall how rejection method works in case q E (—1, 1). 

Applying algorithm described in [10], the rule to get one realization of random 
variable having density fn is as follows. 

(1) we generate two variables: X ^ Je {■) and Y ^ U {0, 1) 

(2) seiT = M{q)fE[X\q)l,fH{X\q). 

(3) If FT > 1 then set Z = X otherwise repeat (1) and (2). 

To see how this algorithm works, we present two simulation results performed 
(consisting of 2000 simulations) with q — —.% (red dots) and q—.% (green dots). 



Fig. 2. Simulation of i.i.d. sequences from P}j. 

Unfortunately this algorithm turns out to be very inefficient for q close to — 1, 
more practically less than say —.85. One can see this by examining Figure 1. Values 
of M (q) are very large then, showing that one needs very large number of obser- 
vations from density fE to obtain one observation from fn- Thus there is still an 
open question to generate efficiently observations from fn for values close to —1. 

One might be inclined to use formula l|2.2p and inversion method applied to its 
finite approximation ad again using procedure 'FindRoot'. Well we applied this 
idea to simulate 2000 observations from Ph for q — —.95. It worked giving the 
following results: 
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Fig3. i.i.d. sequence from Ph by inversion method for q — —.95. 

We used procedure 'FindRoot' of Mathematica. It worked as one can see however 
it lasted quite a time to get the result. 

Besides, when we tried to get 2000 observations from Ph for q = —.97, numerical 
errors seemed to play an important role as one can notice judging from black dots 
that appeared between levels and —.5 on the picture below. 



Fig. 4. i.i.d. sequence from Ph by inversion method for q = —.97. 

4. Proofs 

Proof of TheoremUi Let us denote z = . Hence \z\ < 1. We have 



k=0 

Now let us notice that 

,2 .9 1- i. , l./9\ /, t ^ k/2 



(1 + q'Y - 4zV = (l + q' - ^zq'^'^) (l + + 2^9* 

Now notice that since \z\ < 1, we see that — z'^ + iz\ — 1. Thus we can write 
\/l — z'^ + iz — exp {i9) where 9 — arcsin z and also Vl — = cos 9. Hence we can 
write 

(1 + g'^y - Az^q" = (1 + e^^'q'') (l + e-^^'q'^) , 

and consequently 



^ ' ' 47rcos6l ^ ' ^ 

We will now use so called 'triple product identity' (see |T], Theorem 10.4.1., p. 497) 
that states in our setting, that 



fe(fe-l)/2 2ike 

i 

k— — oo k— — oo 

)/2^-2ije 

k=l J=l 



k=l k=2 
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Now notice that 

(l + e-2«e) ^ 2 cos 61 



cos 9 e*^ cos I 

Hence, 



g-2j(fe-i)e 



\ fc=l fc=2 / 

27r 

(OO CXD N 

1 + e^'' + ^ qi'i) (cos 2fc6l + cos 2 (fc - 1) 61) + i ^ g(') (sin 2/^6* - sin 2 (fc - 1) 6i) 
fc=2 fc=2 > 



COS ( 



^^VgQ) cos(2/fc-l)0. 



fe = 2 

To return to variable z we have to recaU definition of Chebyshev polynomials. 
Namely, we have 

cos ((2fc —1)6')= cos {{2k — 1) arcsin z) 



cos \ {2k — 1) y — — arccosz 
l|_-|^-jfc+i _ 2^ arccos z) 



{-lt+' ^l-Z^U2k-2 {Z); 



where Un{z) is the Chebyshev polynomial of the second kind. More precisely we 
have here: 



(4.1) Un{z) = sin {{n + 1) arccosz) / v 1 — z^. 

It is well known, that sequence {[/„} satisfies three-term recurrence equation 

[/„+l(z) - 2xUn{z) + Un-l{z) = 0, 

with U-i{z) ~ 0, C/q [z) — 1, and can be calculated directly (see Theorem 7.2, 
p. 188) as in (|1.6p . Thus we have shown that 



fH {Z) = ^VT^ (^E(-l)'"'?^^^t/2.-2 {Z)j , 

or equivalently 

fH = ^^v/4-(l-g)x2 (^£(-l)'=-i<z(^)c/2fe-2 (xv/r^/2) j 



□ 



Proof of Corollary\^ We have -Fff (y) = /-2/VT^-^^ "^^^ change vari- 

able to z = — q/2 and use (|2.ip . Thus 

i^/^ (2/) = f Vl^dz+Y.U=2 (-1)'"' 9^^) VT^f/2.-2 {Z) dz). 

We use now the following, easy to prove, formulae: 

VT^U2n {z) d z^y/T ^ {U^n+i {v) / (4n + 4) - U2n~i {v) / (4n)) and /^^ ^T^rfx 
■| + ^ arcsin y + ^y\/l ~ and after rearranging terms get (|2.2p . □ 
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Proof of Corollary ??. Let fu {x\q) = ^^^4- (1 - q)x^ for |a;| y/T~~q < 2. As- 
sertion i) is obtained directly after noting that 

oo 

fH {x\q) = fu{x\q) n ((1 + - (1 - 9) ^'9') (1 - 9') • 
k=l 

Following (|2.ip . we get 

/// (:^;k) = /a (a;k) ( f;(-l)'=-^g(^)c/2fc-2 ) . 



vfe=i 



ii) and mz) are obtained by inserting x — and a; = 1/^1 — q in (|2.1|) and canceling 
out common factors. From (|2.ip it follows also that values U2n (0) and U2n (1/2) 
will be needed. Keeping in mind ()4.ip we see that [/2ti(0) = cos (mr) = (—1)" and 

{1 if n = 3m 
if n = 3m+l . 
— 1 if n = 3m + 2 

On the other hand we see that H^i ((1 + ?'')^ - (1 - '7)a;^g'') n^o(l ~9''^^)L=o 



n (1 - 9^'^) (1 + q") and nr=i ((1 + - (1 - '?)^''z'^) nr=o(l - 



fe=l 



= nr=i((i+'z')'-9')nr=i(i-'z') = nr=i(i+9'+'z'')nr=i(i-9') - 
nr=i ^ nr=i(i - 1') = uZi (i - 1'') ■ 

iv) Putting X = ±-7== in (|2.3p we get 



00 



fe=l A;=l 

Now recall that U2k (1) = 2fc + 1. 

u) To see this notice that g— Hermite polynomials are orthogonal with respect 
to the measure with density fH- Thus we have 

Vm > : / Hm {x\q) fn (x, q) dx = 0. 

J-2/VT^ 

Using (PTTjl know that Vto > 

Lt/cS^™ (^k) V4-(i-g)^^ (1 + Er=2(-i)'+'9(')^2.-2 (xyr^/2)) 

0. Observing that function //f is symmetric and q— Hermite polynomials of odd or- 
der are odd functions, we deduce that above mentioned identities are trivial for odd 



dx 



m. Thus, let us concentrate on even m. Introducing new variable z = x^l — q/2, 
and multiplying both sides of this identity by (1 — 5)™^^ we get 



z^dz = 0, 



(4.2) Vm > 1 : J ^h^ {z\q) (^l + £(-l)'+'g(^) C/2fe-2 (^)j VT 

where hm iz\q) = (1 — g)™^^ Hm {2z/y/l — q\q) . Polynomials hm are called contin- 
uous g— Hermite polynomials. It can be easily verified (following (jl.Sp ) that they 
satisfy the following three-term recurrence equation 

hn+i it\q) = 2thn {t\q) - (1 - q")hn-i {t\q) , 
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with h-i ~ 0, ho {z\q) — 1. Moreover, it is also known that (see e.g. [T]): 



(4.3) hn{cose\q)^Y. 



cos((n- 2k) 9). 



k=0 

Let us change once more variables in (14. 2p and put z = cos r. Then, for Vm > 1 : 
C h^{cosT\q)-K (l + Er=2(-l)''+'9^'^'^2fe-2(cosr)) xsinrdr^O, or 



/ OO \ 

Vm > 1 : y /i™(cosT|g) f ^(-l)'''"ig(2)sin(2fc-l)Tsinrj dr = 
Keeping in mind that 2 sin (2fc — 1) t sin r — cos (2A: — 2) t — cos {2k) t, we see that 
Vm > 1 : ^ /i2m (cosT|g) ^X^'^-l)'""^?'^^^ (cos (2fc - 2) r - cos (2fc) t)^ dr = 0. 
Now keeping in mind that JJ^ cos 2mTdT — for m > we see that 

: ( h2m{cosT\q)\ V (-l)'=+^g(^)(cos(2fc-2)T-cos(2fc)T) )dr = 0. 



Vm > 1 



On the other hand taking into account (|4.3p we see that 

2m 
m — k 



/•TT r 
/ /i2m (cos r|g) COS 2fcr(ir = TT 



for fc = 0, 1, . . . , m, . Hence we have <{2.b\ . 

vi) Keeping in mind that fn is symmetric with respect to x we deduce that 
the point of change of modality of fn must be characterized by the condition 
Ih (Oko) = 0. Calculating second derivative of the right hand side of (|2.ip and 



remembering that I ■\/4 — (1 — (/) 



= we end up with and equation = 

:r— 

- V Er=o -7'=('^+i'/2 + 2(1 -q)x i-'^t Kk + l)(;fe(fe+i)/2(_i)fc+i. Now since 

4:k{k + 1) + 1 = {2k + 1)^ we get equation in vi) defining go- O 

To prove Lemma [T] we need the following lemma. 
Lemma 3. Suppose < r < 1 and n > 3. Then 

T.(^k-l)ri^)<——^. 

k>n ^ ' 

Proof. Recall that for \p\ < 1 we have: X]i>i*P*^^ = (i^) ^ ('^-pf ^^sX 
TZ, P' = '-f^- Th^s we have Y.k>ni'^k l)r(^) = 

^(n-l)(n-4)/2 (2fc - 1) 7-'"('""l)/2-("-l)("-4)/2 = ^(n-l)(n-4)/2 ^ 

Efe>„ (2fc- l)r2"^-2^''-(fe-i)/2-(n-i)(«-4)/2-2(fe-i)^ ^ow uoticc that fc (fc - 1) /2 - 
(71 - 1) (n - 4) /2 - 2(fc - 1) = \ {k- n){k + n~h) > for /c > n > 3. Hence 
Efc>„(2fc-l)r(^) < r("-i)("-4)/2 J2k>n {2k - 1) r2fe-2 < 2r("-i)("-4)/2 ^^^^^^ k (r^) 

o^fa-l)f»-4)/2 "('-')"''-("-l)(''')" < o «r"'"-^)/^ n 



fe-1 
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Proof of Lemma[^ i) 
We have 



Un{cos0) = 



sin((n+ 1)61) 



From this fact we deduce that sup |J7„ (a;)| — {n + 1) . Now using Lemma we 

kl<i 

have 



ri-l 



fe=2 



= V4 - (1 - g) 



< 



F„{y)-l~^ arcsin (i^) - ^-i ,(^)c;^_^ ( 



since sup|^|<]^ Vl — a;^ | C/„ (a::)| = 1. Now to get m) we use routine transformations 
and sum two geometric series. □ 

Proof of Lemma\^ Notice that comparing definitions of fn and fE we have fu {x\q) 

- fE {x\q) X i^i^ X nr=4 ((1 + - (1 - 9)-^?'=) X nr=i (i - - (i + 

g)(l - g3) (1 _ g5) (1 _ ^9) X ^ 1^(1 + ^/c)2 _ _ ^^^2^fej (1 _ ^fe) ^ Now if q e 
[0, 1) we have 

((1 + q'^f - (1 - q)x^q'^) (l - < (l + (l " 

= (l-g2/c) {1 + q'^). 

liqe (-1,0) then 

((l + g'^)'-(l-g)xV) (l-'?'^)< 
, r ((l + g'=)'-4g'=) z/ fcisodd 

^{l-q') {l + \q\')\ 



if k is even 

2 



since then sup 

kl<2/Vl-9 



((1 + q'^f - (1 - g)a;V) - ((1 + - 4g'=) = (1 - 



□ 



4.1. Appendix. Program in Mathematica that generates i.i.d. sequences from //f . 

QA^/^g., M_]:=(Label[pocz]; Y=y/. FindRoot[F[y, q]-RandomReal[], { y, 0} ]; 

u=RandomReal[];t=newMM[q,M]/R[Y, q,Mj;If[t u< =1, Y, Goto[pocz]]); 

However it requires definition of function F which is in fact function Fe of this 
paper. It is quite lengthy. newMM denotes function M of this paper. Further 
function R denotes the ratio fn/fE- Parameter M denotes number that we insert 
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instead of oo in the above mentioned formulae. The above procedure produces 1 
observation from fn- 

Now AA[q_,M_,h_J--^ListPlot[Table[QN[qM],{2000}],PlotStyle->Hue[h]]; produces 
table of 2000 observation from fni-lq) and plots it in color h. Then AA[.8,100,.4] 
and AA[-.8,100,0] produce plots for g = .8 in color. 4(green) and g = — .8 in color 
O(red). 
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